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Asymmetric nuclear matter in a parity doublet model with hidden local symmetry 
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We construct a model to describe dense hadronic matter at zero and finite temperature, based 
on the parity doublet model of DeTar and Kunihiro, with including the iso-singlet scalar meson a 
as well as p and ut mesons. We show that, by including a six-point interaction of a meson, the 
model reasonably reproduces the properties of the normal nuclear matter with the chiral invariant 
nucleon mass mo in the range from 500 MeV to 900 MeV. Furthermore, we study the phase diagram 
based on the model, which shows that the value of the chiral condensate drops at the liquid-gas 
phase transition point and at the chiral phase transition point. We also study asymmetric nuclear 
matter and find that the first order phase transition for the liquid-gas phase transition disappears in 
asymmetric matter and that the critical density for the chiral phase transition at non-zero density 
becomes smaller for larger asymmetry. 
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I. INTRODUCTION 

With the advent of next generation radioactive beam 
facilities isospin asymmetric nuclear matter claims much 
attention in contemporary nuclear physics. At those fa¬ 
cilities we could create terrestrial environment to study 
dense matter with a large neutron or proton excess 
through nuclear reactions with radioactive nuclei. 

Studying nuclear matter is also important to under¬ 
stand the structure of neutron stars [1]. In 2010 and 2013, 
two neutron stars with twice solar mass were found [2, 3] 
and many models yielding the soft equation of states 
(EOS) were excluded. Neutron stars offer very cold and 
asymmetric dense environment and may have hyperons 
in the core of the stars. If there are hyperonic degrees 
of freedom, it is expected that the EOS becomes softer 
and neutron star mass becomes lighter. Another impor¬ 
tant astrophysical site for nuclear matter is a hybrid star 
whose center has quark matter [4]. 

The properties of asymmetric matter have been in¬ 
vestigated in various approaches [5-14]. Very recently 
liquid-gas and chiral phase transition are studied in a par¬ 
ity doublet model with a six-point scalar interaction in 
which mesonic fluctuations are included by means of the 
functional renormalization group [15]. In this work we 
study isospin asymmetric dense matter in the framework 
of the parity doublet model (mirror assignment) [16, 17]. 
The properties of symmetric dense matter such as chiral 
phase transition were extensively studied in the parity 
doublet models at zero or finite temperature [18-22]. We 
extend the parity doublet model by including p and w 
mesons through the hidden local symmetry and also by 
adding a six-point interaction of a scalar meson. Here, 
as a first step, we will not consider hyperonic matter and 
work within the mean field approximation. 

We determine our model parameters, except the chi¬ 
ral invariant mass (mg), by performing global fitting to 
physical inputs (masses and pion decay constant in free 
space and nuclear matter properties). We then study the 
equation of state and the phase diagram of dense mat¬ 


ter at finite temperature. We find that the predicted 
slope parameter at the saturation density meets the con¬ 
straint from heavy ion experiments and neutron star ob¬ 
servations (see, e.g. Refs. [23, 24]) and observe that the 
chiral condensate drops at the chiral and liquid-gas tran¬ 
sition points. It is also seen that smaller mg values prefer 
smaller critical densities for chiral phase transition. The 
study of asymmetric matter reveals that the first order 
nature of the liquid-gas transition disappears in asym¬ 
metric matter and the critical densities for the chiral 
transition become smaller with increasing asymmetries, 
which are consistent with previous studies. 

In section II we extend the parity doublet model, and 
in section III we fix the model parameters. Our results on 
bulk properties of nuclear matter and density dependence 
of chiral condensate and nucleon mass are given in section 
IV. We present the phase diagram of dense (asymmetric) 
matter in section V. Finally, conclusion and discussion 
follow in section VI 


II. EXTENDED PARITY DOUBLET MODEL 

We construct a chiral effective model based on the par¬ 
ity doublet model [16, 17], in which a nucleon with posi¬ 
tive parity is regarded as a chiral partner to the one with 
negative parity and they belong to the same multiplet. 
The transformation properties of the positive and nega¬ 
tive parity nucleon fields are given by 

V'lr -t Ipll -t OL'tpll (2.1) 

V'2r gL'02r, t/’2/ gR'tl^2l , (2.2) 

where qr {pl) is an element of SU(2)fl (SU(2 )l) chiral 
symmetry group, and ifir and '!/' 2 r {i’ll and 'ip 2 i) are the 
right-handed (left-handed) fields projected as 

1plr,2r = Pr'4’1,2 , i’ll,21 = Pl ' 01,2 , 
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(2.3) 



2 


To construct a linear sigma model including these nu¬ 
cleon fields, we introduce the following field: 

M = a + in ■ f , (2-4) 

where a denotes an iso-singlet scalar field, tt the pion 
field, and r the Pauli matrices. The chiral transformation 
property of M is given by 

M ^ gLMg^^ . (2.5) 

By using these fields, the Lagrangian of a linear sigma 
model type is expressed as ^ 

Cn ='4’irh^Dfj,ipir -b 

+ i’2ril^D^'lp2r + '4>2lil^ D 

- Too [tplli’2r - i’lri’21 “ ^2li’lr + i’2ri’ll] 

- gi + -ipiiM-tpir] 

- 92[^2TM%l^2l+4’2lM'^i’2r] , (2-6) 

where mg is the chiral invariant mass, gi and g 2 are the 
coupling constants, and the covariant derivatives include 
the external gauge fields TZ^ and as 

D^i^lr,2l = {dfi — iTZfj_) V'lr,2/ , 

2r — V’l/,2r ■ (2T) 

The meson part of the Lagrangian is given by 

Cm = Jtr [d^Md^M^] -Vm - Vsb , (2.8) 

where the meson potential Vm and the explicit chiral 
symmetry breaking potential Vsb are 

Vm=- ^/f^tr [MMt] -b {tr [MM^] f 

-^A6{tr[MMt]}^ (2.9) 

f/sB=-ie(tr[A4l'M] -btr[A4Ml']) . (2.10) 

Here A4 is the quark mass matrix given as 

and £ is a constant of mass dimension two. In the present 
analysis, we neglect the isospin breaking effect due to 
the mass difference of up and down quarks and take 


^ This Lagrangian is rewritten into a more familiar form in the 
literatures as 

+ '02*7''d^,i/’2 - rno (■0175'i/’2 - '0275V’l) 

— gi'tfii (it + ijsn ■ f) ifii — g 2'02 (t — 1757 ? • r) -(/;2 . 


niu = rud = m. It should be noticed that the poten¬ 
tial Vm includes the dimension-six term which will play 
a very important role in reproducing the properties of the 
normal nuclear matter with rather wide range of the chi¬ 
ral invariant nucleon masses. (See next section.) When 
Ae > 0, the potential Vm + Vsb is not bounded from be¬ 
low. In the present analysis, we determine the vacuum 
in the following way: We first solve the stationary con¬ 
dition of the potential. When there are more than one 
solutions, we choose the one with the lowest energy at 
the stationary point. 


We next include p and to mesons into the model based 
on the hidden local symmetry (HLS) theory [25, 27]. The 
HLS is introduced by performing the polar decomposition 
of the field M as 


M = (2.12) 

where a is a scalar meson field and ^nd transform 
as 


^L,R -t hu:hp^L,RgL^R (2-13) 

with huj € U(I)jjBs ^ ^U(2 )jjls- In. the unitary 

gauge, and are parameterized as 

^R=^l= expizn-TyU) , (2.14) 

where T® = ra/2 (a = 1,2,3) with being the Pauli 
matrix. In the HLS, the vector mesons are introduced as 
the gauge bosons of the HLS which transform as 

cjp y • ojp • -b • /fl, (2.15) 

Pp ^ hp ■ pp ■ hl +—dphp-, (2.16) 

9p 

where g^^ and gp are the corresponding gauge coupling 
constants. 


To construct a model Lagrangian with the HLS, it is 
convenient to introduce the following one-forms: 




(2.17) 


where the covariant derivatives are given as 


= d>^^L + igpP^^L + igojOO^^L + i^LC^ , 

D^^r = + igpP^^R + igujijJ^^R + i^rTZ^ ■ (2.18) 


Now, the mesonic part of the Lagrangian extended by 
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where the first line is from Eq. (2.8) with the following 
form of the potential: 

Va = —— 0^6(7® (2.20) 

VsB = — ^'hiecr tr [U + C/^] . (2-21) 

The second and third lines contain the mass and kinetic 
terms of vector mesons, respectively. 

We also rewrite the nucleon part of the Lagrangian as 

Cn + 'lp2Til^ D ^'lp2T + 

- Too [V'l/V’2r - 4>lr'tp2l “ 4>2l^lr + i’2ri’ll] 

- 9lO- + TpllU^ir] - 920- ['tp2rU'tp2l + tp2lU'‘1p2r] 

- apNN ■0ii7^(?ia||,iCL)V'i/ + i’lrl^{^Ra\\f,^R)'tplT 

- QpNN '02i7^(C/jQ:||M?fl)^2; + ^2r7^(?ia||/iCL)V'2r 

- aoNN tr [dll^] + ■iplrl^lplr + 'tp2a^^2l + i’2rl^i’2r) 


( 2 . 22 ) 


the HLS is expressed as 

Cm =^dpad^a + cr^tr [a±pa'^] - - Vsb 



(2.19) 


The vacuum expectation value (VEV) of the cr field, 
denoted by (Jq, is determined by the stationary condition 
for the potential, as we explained above. The non-zero (Tq 
breaks the chiral symmetry spontaneously and generates 
the masses of nucleons as 

r,„„. ^ CO ^ (2.23, 

We obtain the masses of the positive-parity and negative- 
parity nucleons by diagonalizing the mass matrix. Here, 
we write the mass eigenstates as N+ and N-, which are 
related to and ■i /'2 as 

fN+\-f cos 6 » 75sin6»\ /V'lO 

\N_) ~ ^^-yssind cos 6 ) ^^'27 ’ ^ ’ 

where 9 is the mixing angle given by 

tan 20 = -—— . (2.25) 

(ffi + 02)0-0 

The mass eigenvalues are determined as 

w± = ^ (^ 1/(51 + 02 ) -k4m§ T (01 - 02 ) 0 - 0 ^ , 

(2.26) 

where m+ and m_ are the masses of positive and nega¬ 


tive parity baryons, respectively. ^ From this expression, 
one can easily see that the spontaneous chiral symme¬ 
try breaking is responsible for the mass differences of the 
parity partners. 


III. DETERMINATION OF MODEL 
PARAMETERS 

In this section, we determine the 10 unknown parame¬ 
ters in this model by performing a global fit with chosen 
mo to masses and the pion decay constant in free space, 
and to normal nuclear matter properties. In this global 
fitting we have used mo = 900,800, 700,600 and 500 MeV 
as we have found that with mo = 400 MeV or less we can¬ 
not reproduce the normal nuclear matter properties such 
as incompressibility. 

Determined parameters are summarized in Table I. 
Now, we describe how we fix the model parameters and 
what are the inputs. 

First, we determine 6 parameters by using the phys¬ 
ical inputs listed in Table II. We choose the mass of 
positive- (negative-) parity nucleons as m+ = 939 MeV 


^ As we will explain in the next section, we use = 939 MeV 
and m_ = 1535 MeV as inputs to determine the values of 31 and 
32 for given mo. When we take as the mass of the negative 
parity baryon and m_ as the one of positive parity baryon, i.e. 
m_ = 939 MeV and m+ = 1535 MeV, the determined values of 
31 and 32 are exchanged. One can verify that this exchange does 
not cause any physical difference by swapping ipi with r/ 12 . 
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TABLE I: Determined model parameters for given mo. Here 
rriui = 783 MeV, rUp = 776 MeV and me = 


mo [MeV] 

500 

600 

700 

800 

900 

5i 

15.4 

14.8 

14.2 

13.3 

12.3 

52 

8.96 

8.43 

7.76 

6.94 

5.92 

OuiNN 

11.4 

9.12 

7.31 

5.67 

3.54 

PpNN 

8.05 

6.97 

7.46 

7.75 

8.75 

/7[MeV] 

435 

434 

402 

316 

109 

A 

40.5 

39.4 

34.5 

22.5 

4.26 

-^6 

16.3 

15.4 

13.5 

8.66 

0.607 


where the asymmetry parameter S is defined as 

<5=^^^ = ^. (3.5) 

Pb Pb 

We fit the remaining four-parameters to the empirical 
values shown in Table III. 


TABLE III: Physical inputs. = 923MeV. We took the 
value of incompressibility and symmetry energy from [28] and 
[23] respectively. 


po{p*B)[im ^] 

F;/A(^*fj) -m+[MeV] 

A [MeV] 

Esym[MeY] 

0.16 

-16 

240 

31 


(m_ = 1535 MeV). As in some literatures, one may also 
try m_ = 1200 MeV, but we have chosen the lightest 
and observed one. The pion mass and decay constant 
are = 140 MeV and = 93 MeV. 


TABLE II: Physical inputs in vacuum (MeV). 


m+ 


rricp 

rUp 

U 


939 

1535 

783 

776 

93 

140 


We would like to stress that we can reproduce the value 
of the incompressibility by the inclusion of the six-point 
interaction of the scalar meson ct, in contrast to the pre¬ 
vious analyses by parity doublet models [19, 20] where it 
seems difficult to reproduce the small value of the incom¬ 
pressibility. 

IV. EQUATION OF STATE 


Next, the remaining parameters are fixed by the satu¬ 
ration density, the binding energy, the incompressibility 
and the symmetry energy for the normal nuclear matter 
at zero temperature. The normal nuclear matter proper¬ 
ties for the fit are discussed below, and empirical values 
of them are summarized in Table III. 

The typical value of the nuclear saturation density and 
binding energy are 0.16 fm“^ and —16 MeV, respectively. 
With TO+ = 939 MeV we have 


p{p*B = 923 MeV) = po = 0.16 fm-^ 


'E 


€ 

A ■ 


-m+ 

Po 

.PB 


-16 MeV . 


(3.1) 

(3.2) 


We can consider another property of nuclear matter, the 
incompressibility, which is given by the curvature of bind¬ 
ing energy at saturation density and corresponds to the 
“hardness” of the matter: 




dp^ 


= 9po 


dpB 


Po 


dp 


(3.3) 


Po 


The symmetry energy per nucleon is from the difference 
of proton and neutron and given as 


^sym{PB^ — 


1 d^iE/A) 
2 ! dS^ 

1 dHe/p) 

2 ! 9(52 ’ 


(3.4) 


In this section, we study the equation of state for cold 
nuclear matter using the model constructed in the previ¬ 
ous sections. 

In Fig. 1, we show the dependence of the binding en¬ 
ergy (upper panel) and the pressure (lower panel) on the 
baryon number density for mo = 500 MeV as an exam¬ 
ple. The red line in the upper panel of Fig. 1 shows the 
dependence of the binding energy on the baryon number 
density of symmetric matter. From this, one can eas¬ 
ily see that the binding energy is actually minimized at 
Pb = Po, which implies that the existence of the bound 
nuclei and the saturation property are actually repro¬ 
duced.^ 

The lower panel shows that the pressure for Pn/ Pb = 
0.5 is negative in the low density region, which means 
that the liquid phase of hadron coexists with the gas 
phase. In asymmetric matter, on the other hand, when 
the degree of asymmetry, Pn/ Pb, is around 0.6~0.7, the 
coexistence phase disappears and the pressure increases 
monotonically for larger asymmetry, which implies that 
the bound nuclei do not exist for Pn/ Pb > 0.7. In 
case of the pure neutron matter {pn/P b = 1 )) the energy 
density as well as the pressure monotonically increases 
with the baryon number density. The EoS (pressure) of 


® We would like to stress that we do not use the minimization 
condition as an input. 

^ Actually, the critical value of Pn/pB depends on the chiral in¬ 
variant mass: pn/P b ~ 0.7 for mo = 500 MeV and pn/P b ~ 0.8 
for mo = 900 MeV. 
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FIG. 1: Density dependence of the binding energy (upper 
panel) and the pressure (lower panel) for mo = 500 MeV. ps 
is the baryon number density and is the neutron number 
density: pn/pB = 0.5 implies symmetric nuclear matter and 
pn/pB = 1 pure neutron matter. 


asymmetric nuclear matter is also discussed in molecular 
dynamics [29] and many body perturbation [30]. Though 
our analysis is done within the mean field approximation, 
the qualitative behaviors in symmetric and neutron mat¬ 
ter from [29, 30] are similar to our result. 

Next, we study the slope parameter, L, which is the 
gradient of the symmetry energy at the saturation den¬ 
sity: 


L — 3po 


d^symiPB') 


dpB 


(4.1) 


The constraint for the value of L is obtained from some 
experiments, for example, heavy ion collision experi¬ 
ments, nuclear masses and so on [23, 24]. In Table IV, 
we list our calculated values of the slope parameters for 
several choices of the chiral invariant mass niQ. This 
shows that the slope parameter hardly depend on mo, 
all of which are within the allowed region shown in 
Refs. [23, 24]. 

In Fig. 2, we show the dependence of the VEV of the 
a held (upper panel) and the baryon number density ps 
(lower panel) on the baryon number chemical potential 
Pb for Too = 500 MeV in symmetric matter (red curves) 
and in the pure neutron matter (green curves). 


TABLE IV: Predicted values of the slope parameter for mo = 
500 MeV to 900 MeV. 


mo [MeV] 

L[MeV] 

900 

75 

800 

74 

700 

78 

600 

78 

500 

75 



fielMeV] 



^BlMeV] 


FIG. 2: Density dependence of (tq (upper panel) and chem¬ 
ical potential dependence of ps (lower panel) for mo = 500 
MeV and pi=0. 


The upper panel shows that there are two points, 
Pb 900 MeV and 1300 MeV, where the value of cto 
changes rapidly. In the case of symmetric matter shown 
by the red curve, the hrst jump for pB ^ 900 MeV can 
be identihed with the hrst order phase transition from 
the liquid phase to the gas phase as the baryon number 
density (an order parameter of the liquid-gas transition) 
undergoes sudden change around pb ^ 900 MeV (see 
the lower panel of Fig. 2). ® Then, the phase transition 


® One may schematically understand why the chiral condensate 
drops at the liquid-phase transition, where the baryon num¬ 
ber density jumps as we increase the baryon chemical potential, 
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around ^ 1300 MeV can be naturally identified as 
the chiral phase transition. In symmetric nuclear mat¬ 
ter (red dashed-line), the liquid-gas phase transition is 
first-order and there exists coexistence phase. The exis¬ 
tence of coexistence phase in the nuclear liquid-gas phase 
transition has been confirmed in the experiments [31]. 
The coexistence phase in symmetric matter disappears 
by increasing /r/, and the liquid-gas transition becomes 
second order as suggested in [6]. 





FIG. 3: Density dependence of the effective nucleon masses 
for mo = 500 MeV at p/ = 0. 

Figure 3 shows the density dependence of the effective 
mass of positive (negative) parity nucleon, m+ (m_). As 
Pb increases, m+ and m_ gradually get close to the chi¬ 
ral invariant mass. This is a feature from parity dou¬ 
blet structure, and two nucleon masses degenerate to toq 
when the chiral symmetry is completely restored. How¬ 
ever, the mass difference between the parity partners re¬ 
mains finite in our model due to the current quark mass. 
The critical density for chiral symmetry restoration de¬ 
pends also on the chiral invariant mass, which will be 
discussed in detail in the next section. 


V. PHASE DIAGRAM 

In this section, we explore the phase structure of our 
model at finite temperature and density with the isospin 
asymmetry. Our primary interest here is to see how the 
onset of the chiral and liquid-gas phase transition de¬ 
pends on the isospin asymmetry and the chiral invariant 


through the Pauli exclusion principle. As the baryon density in¬ 
creases, the low-lying phase space relevant for quark-antiquark 
condensates is occupied by the fermions (quarks in nucleons in 
this case) that constitute the Fermi sea; therefore, forming quark- 
antiquark condensates requires much energy. At the liquid-gas 
transition point, there is a sudden increase in the number den¬ 
sity, and so we could expect that the chiral (quark-antiquark) 
condensate changes drastically. 


nucleon mass mo- In this study, we will not consider the 
charged pion condensation, and so we take |/i/| < 



fielMeV] 



tielMeV] 


FIG. 4: Phase diagrams for mo = 900 MeV at p/ = 0 
(upper panel) and pi — 100MeV (lower panel). The solid 
line is for the first order phase transition, the dashed line for 
the crossover, and the point for the critical point (2nd order). 

The upper panel of Figure 4 corresponds to the phase 
diagram ai pi = 0, where we have the first order 
liquid-gas phase transition (blue solid line) around pb = 
900 MeV and chiral phase transition (red solid line) 
around pB = 1500 MeV at zero and finite but small tem¬ 
peratures. The magenta dots are second order critical 
points, and green dashed line shows the crossover. 

The two crossover lines meet at a point around 
{pb, T) = (700 MeV, 100 MeV) as seen in a study with 
a parity doublet model [20]. 

The lower panel of Figure 4 shows the phase diagram 
of asymmetric matter with mo = 900 MeV at pi = 100 
MeV. The trend of the liquid-gas and chiral phase tran¬ 
sition is almost the same with the upper panel {pi = 0); 
however, the liquid-gas transition becomes second order 
even at zero temperature as briefly mentioned in the pre¬ 
vious section. 

Now, we study how the chiral invariant mass affects 
the phase diagram. We take mo = 500 MeV as an exam¬ 
ple. Since the chiral invariant mass measures the amount 
of spontaneous chiral symmetry breaking needed for the 
nucleon masses and their mass splitting, we may expect 
that it mainly affects chiral symmetry in the phase dia¬ 
gram. The results are shown in Fig. 5. As expected, the 
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nature of the liquid-gas transitions does not change, while 
the first order chiral phase transition becomes second or¬ 
der. Furthermore, as mg decreases, the critical density 



tiglMeV] 



fielMeV] 


FIG. 5: The phase diagrams for mo = 500 MeV. 

(chemical potential) for chiral transition decreases mono- 
tonically as summarized in Table V. Note that the sat¬ 
uration density is an input, and so the corresponding 
chemical potential, in our study. This is why the critical 
density (chemical potential) for liquid-gas phase transi¬ 
tion does not change. 


TABLE V: Critical values of the baryon chemical potential 
a.t T — 0 and fii = 0. 


mg [MeV] 

PBfg[MeV] psFjfm 

MBy[MeV] 

PBy[fm-3] 

900 

923 

0.16 

1540 

2.98 

800 

923 

0.16 

1643 

1.75 

700 

923 

0.16 

1554 

1.09 

600 

923 

0.16 

1426 

0.69 

500 

923 

0.16 

1305 

0.44 


When we employ mo = 900 MeV, the critical den¬ 
sity for chiral transition is about 17po- On the other 
hand, the critical density becomes ~ 3po with mg = 500 
MeV. Comparing with a previous study in a parity dou¬ 
blet model with mg = 800 MeV [20], we find that our 
result for the symmetric matter is almost the same with 


the one from the previous study. Note, however, that we 
can explore the phase diagrams with the chiral invari¬ 
ant mass of the range 500 MeV < mg < 900MeV as we 
introduced the six-point interaction. 

Table VI shows the critical baryon chemical potential 
and density for /i/ = 100 MeV; comparing with Table 
V, we observe that the critical values become (slightly) 
smaller. Here the liquid-gas phase transition is second 
order so that psfg must be 0. 


TABLE VI: Critical values of the baryon chemical potential 
and density at T = 0 and pi = 100 MeV. 


TOg [MeV] 

PbU^bV] psrjfm-3] 

MB^[MeV] 

PB^[fni-3] 

900 

891 

0.0 

1537 

2.98 

800 

891 

0.0 

1637 

1.74 

700 

891 

0.0 

1543 

1.09 

600 

888 

0.0 

1410 

0.68 

500 

881 

0.0 

1285 

0.43 


VI. SUMMARY AND DISCUSSION 

We have constructed a model for asymmetric nuclear 
matter by extending the parity doublet model. We in¬ 
troduced vector mesons (p and ui) through hidden local 
symmetry and also included the six-point interaction of 
of a meson. We fixed our model parameters with chosen 
mg by performing a global fit to physical inputs (masses 
and pion decay constant in free space and nuclear matter 
properties). With the six-point potential, we were able 
to reproduce normal nuclear matter properties with mg 
in the range from 500 MeV to 900 MeV. 

We first studied the equation of state and the phase 
diagram of dense symmetric matter at finite tempera¬ 
ture. We observed that the slope parameter at the sat¬ 
uration density satisfies the constraint from heavy ion 
experiments and neutron star observations, for instance, 
see [23, 24] and found that the chiral condensate changes 
drastically at the chiral and liquid-gas transition points. 

We then moved on to asymmetric dense matter with 
non-zero iso-spin chemical potential. We showed that 
the first order nature of the liquid-gas transition disap¬ 
pears in asymmetric matter and the critical densities for 
the chiral transition becomes smaller with increasing iso¬ 
spin chemical potentials, which are in agreement with the 
results from existing literatures. We also showed that 
smaller chiral invariant nucleon mass favors smaller crit¬ 
ical density for chiral phase transition both in symmetric 
and asymmetric dense matter. 

In this study we don’t consider interesting phenomena 
in dense (asymmetric) matter such as transition from nu¬ 
clear matter to hyperonic matter, charged pion conden¬ 
sation with large iso-spin chemical potential, and so on. 





























These will be relegated to future works. 

Finally, we discuss the chiral invariant nucleon mass. 
In our work, we chose the chiral invariant mass mo = 500- 
900 MeV to reproduce the properties of normal nuclear 
matter. However, there are other choices from various 
studies. In [19, 20], nuclear matter was studied in a par¬ 
ity doublet model and rather large value of the chiral 
invariant mass was used mo ~ 800 MeV. On the other 
hand, in [17] the authors determined mo from the decay 
width of N* —> Ntt to be mo = 270 MeV, while in [32] 
they used the decay modes of N* Ntt and Oi —?■ Try 
to obtain mo ~ 500 MeV. Recently, the study of parity 
doublet structure using lattice QCD [33] shows that the 
positive parity nucleon mass changes very small near the 
deconfinement transition, which may imply that mo does 
exist in nature and its value is close to the positive par¬ 
ity nucleon mass. The existence of chiral invariant mass 
is also discussed in the context of Skyrmion crystal [34]. 
Since the values of the chiral invariant mass are diverse in 


literatures, it is quite important and interesting to nar¬ 
row down the value of ttiq and to dig further down to the 
role of mo in hadron physics. 
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